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1.0  INTRODUCTION 

A series  of  optical  homodyne  experiments  (Refs.  1 through  4)  during  the  period 
1962-1964  demonstrated  the  feasibility  of  exploiting  the  Doppler  shift  of  scattered  laser 
light  for  measuring  the  velocity  of  moving  matter.  Of  particular  importance  to  the 
aerodynamicist  was  the  demonstration  by  Yeh  and  Cummins  in  1964  (Ref.  4)  that  fluid 
velocities  could  be  inferred  by  measuring  the  Doppler  shift  frequency  of  laser  light 
scattered  by  particulate  matter  entrained  in  the  fluid. 

Since  1964,  worldwide  interest,  intensive  development  efforts,  and  increasingly 
refined  applications  have  resulted  in  significant  improvements  in  laser  Doppler  velocimeter 
(LDV)  hardware  and  techniques  as  well  as  in  a heightened  awareness  of  remaining 
technological  deficiencies.  In  recent  years,  significant  efforts  have  been  directed  toward 
specifying  and  reducing  the  magnitudes  of  the  systematic  and  random  errors  contributing 
to  the  overall  uncertainty  of  LDV  measurements  (Ref.  5).  One  subtle  systematic  error 
mechanism  not  yet  widely  recognized  is  that  caused  by  not  accounting  for  or  improperly 
accounting  for  axis  nonorthogonality  in  multiaxis  LDV  systems. 

This  error  mechanism  was  recently  encountered  in  a large  transonic  wind  tunnel 
application  of  a two-axis  LDV  system  (Ref.  6).  This  system  was  used  for  free-stream, 
pitch-plane  velocity  magnitude,  and  flow-angle  measurements.  Orthogonal  axis,  data 
reduction  equations  were  initially  used;  however,  the  measurement  axes  were 
subsequently  found  to  subtend  an  angle  of  about  89.1  deg.  This  oversight  resulted  in 
velocity  magnitude  errors  as  large  as  0.8  percent  of  measured  value  and  flow-angle 
measurement  errors  as  large  as  0.4  deg.  The  tunnel  test  requirements  demanded  overall 
uncertainties  in  velocity  magnitude  and  flow-angle  measurements  of  less  than  ±1.0 
percent  of  measured  value  and  less  than  ±0.1  deg.  Thus,  the  oversight  errors  were 
obviously  unacceptable  and  were  ultimately  eliminated  by  using  the  equations  developed 
in  this  report. 

The  measurement  axes  of  virtually  all  two-  and  three-axis  LDV  systems  are  either 
inadvertently  or  intentionally  nonorthogonal  to  some  degree.  Also,  the  LDV  axes  may  be 
rotated  with  respect  to  a reference  set  of  rectangular  Cartesian  coordinates.  Whether  or 
not  special  data  reduction  equations  are  required  to  account  for  axis  nonorthogonality 
depends  upon  the  magnitude  of  the  nonorthogonality  and  the  accuracy  sought  in  the 
particular  LDV  application.  When  special  equations  are  necessary,  it  is  essential  to 
incorporate  into  the  transformation  relationships  the  subtle  point  that  the  LDV  system 
measures  velocity  vector  projections  rather  than  components  when  its  measurement  axes 
are  oblique. 
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At  least  two  previous  analyses  have  dealt  with  nonorthogonal  LDV  measurement 
axes.  The  first  (Ref.  7)  treats  a special  3-space  case  of  transmitter- receiver  symmetry  and 
is  not,  therefore,  universally  applicable.  The  second  (Ref.  8)  provides  only  the  2-space 
transformation  expressions  relating  two-axis  LDV  measurements  to  reference  axis 
components.  A general,  comprehensive,  3-space  analysis  of  LDV  measurement  axis 
nonorthogonality  and  rotation  apparently  is  not  yet  available. 

The  purpose  of  this  report  is  to  derive  the  general  3-space  equations  for  velocity 
vector  magnitude  and  spatial  orientation  in  terms  of  oblique  axis  parameters  and  the 
velocity  vector  projections  measured  by  each  LDV  axis.  This  is  accomplished  by  first 
employing  fundamental  tensor  relationships  to  derive  expressions  that  interrelate  vector 
measure  numbers,  components,  and  projections.  These  expressions  are  then  used  along 
with  the  metric  tensor  to  develop  the  desired  equations. 

2.0  VECTOR  MEASURE  NUMBERS,  COMPONENTS,  AND  PROJECTIONS 

Figure  1 depicts  a three-dimensional  space  simultaneously  covered  by  the  rectangular 
Cartesian  coordinate  system  x*  and  the  general  or  curvilinear  coordinate  system  x'\ 


Figure  1.  3-space  coordinate  systems  x\x2,x3  and  xv,x2',x3 
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Let  x>  and  x’’  be  functionally  related  so  that  the  transformation 

xr  = x’Xx1).  •'  = 1'.  2',  3',  and  i = 1,  2,  3 
is  single  valued  and  reversible.  Let  the  inverse  transformation 

x * = xKx*'),  i = 1,  2,  3,  and  i'  = 1',  2',  3' 


(1) 


(2) 


also  be  single  valued.  Under  these  conditions  the  variable  point  P shown  in  Fig.  1 can  be 
uniquely  specified  in  terms  of  either  its  x‘  or  x’'  coordinates. 

— ► 

Let  r denote  the  position  vector  of  the  point  P with  respect  to  0,  the  origin  of  the 
rectangular  Cartesian  coordinate  system.  A differential  change  in  r can  be  expressed  by 


d?  = , dx* 

(V 


(3) 


where  the  Einstein  summation  convention  is  used.*  For  unit  displacement  along  the  x*' 
axis,  the  change  in  r is  tangential  to  x’'  and  is  equal  to  dr/9x‘\  The  vectors 


dr 

dx *' 


(4) 


are  the  unitary  vectors  with  respect  to  the  point  P (Ref.  9)  which,  in  general,  constitute  a 
nonorthonormal  basis  for  all  field  vectors  defined  at  the  point  P. 

The  unitary  set  defines  a parallelepiped  whose  volume  is  given  by  the  scalar  triple 
product  of  these  vectors  (Ref.  9) 


0 = tr 


*3 ' 


(5) 


’Summation  over  the  index  for  all  admissible  values  is  to  be  carried  out  whenever  the  same  index  occurs 
both  as  a superscript  and  as  a subscript  in  a single  term.  A superscript  in  the  denominator  of  a term  is 
regarded  as  a subscript  on  the  quotient.  The  present  analysis  is  limited  to  3-space  so  that  the  admissible 
values  are  always,  1 , 2,  and  3. 
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A second,  generally  nonorthonormal,  set  of  basis  vectors  can  be  defined  in  terms  of  the 
unitary  set  and  their  scalar  triple  product  (Ref.  9) 

£(?,-*•)  (6) 

These  are  the  reciprocal  unitary  vectors.  The  unitary  and  reciprocal  unitary  vectors  are 
seen  to  satisfy  the  condition 


• V = 


where  d\,  is  the  Kronecker  delta.* 


The  differential  change  in  the  vector  7 can  now  be  expressed  either  in  terms  of  the 
unitary  basis  as 

<8> 

or  in  terms  of  the  reciprocal  unitary  basis  as 

2=  ti'dxr  <9> 

Analogously,  any  field  vector  V defined  at  the  point  P can  be  resolved  into  components 
either  with  respect  to  the  unitary  basis  as 

V = aj ' v 1 (10) 

or  with  respect  to  the  reciprocal  unitary  basis  as 

v = 3*v  (in 

where  v’’  are  the  contravariant  measure  numbers  of  V and  Vj'  are  the  covariant  measure 
numbers.  Equations  (10)  and  (11)  are  statements  of  the  familiar  parallelogram  law  for 
formation  of  the  vector  V by  vector  addition  of  its  components. 

The  squared  length  of  dr  is  expressed  with  respect  to  the  unitary  set  by 

(ds')2  ardxj  dxi'  = g,r  dx‘  dxi'  (12) 


= !,«'=  i' 


s\'  = 0,  i'  4 V 
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where  gj'j'  is  a symmetric  tensor  termed  the  covariant  metric  tensor  (Ref.  10).  This  same 
length  is  expressed  with  respect  to  the  reciprocal  unitary  set  by 


(ds')2  = dr  • dr  = a'  • a(  dxj  - dx^  * = g‘  1 dx  j'dxj ' 


(13) 


where  g1'*'  is  a symmetric  tensor  termed  the  contravariant  meiric  tensor.  The  analogous 
expressions  for  a general  field  vector  V are 


V2  = o . /.  /V1  vl 


giVv  v 


(14) 


and 


V2  = g1  ) vj/vj- 


(15) 


A set  of  unit  vectors  collinear  with  the  unitary  set  is  given  by 
u.  ^ = 1 (no  summation  on  i') 


VfTT. 


(16) 


The  magnitudes  of  the  components  of  V with  respect  to  this  unit  set  are  given  by 


(V).,  = v1  ai ' • 


v gjr  v< 


~ (17) 

— ► 

The  orthogonal  projections  of  V onto  the  axes  defined  by  the  unit  set  are  given  by 


(Vi)j'  = Vj'a^ 


. ' vJ 


£.  ' V 

J 


V si'i ' V «i  'i  ■ 


(no  summation  on  i ') 

v «iV  (18) 


Study  of  Eqs.  (5)  and  (6)  shows  that  the  aj'  and  a)  are  identical  quantities  when  the 
cij'  are  the  orthonormal  basis  vectors  of  a rectangular  Cartesian  coordinate  system. 
Under  these  conditions 


giy  - *.r  - *‘r-  s'r 


(19) 


and  consequently 


v = vi  = (V),*  = (V,v 


CO) 


However,  Vj'f  v*  , (V)j',  and  (Vx)/  are  distinctly  different  quantities  when  the  a*'  are 
nonorthonormal.  It  will  be  shown  subsequently  that  an  LDV  system  measures  the 
quantity  (V^jp. 
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In  Riemannian  geometric  space  the  element  of  arc  length  ds*  is  given  by  Eqs.  (12 
and  13).  If  there  exists  a coordinate  transformation 


(x‘) 


(21) 


such  that 

(ds')2  = (ds)^  = dx1  dx^  (22) 

then  the  Riemannian  space  is  said  to  be  Euclidean  (Ref.  11).  In  Euclidean  space,  any 
coordinate  system  for  which  the  elements  of  the  metric  tensor  are  constants  is  termed  a 
Cartesian  coordinate  system  (Ref.  11). 

All  essential  metric  (quantitative  measures  of  length,  angles,  etc.)  properties  of 
Euclidean  space  can  be  completely  specified  by  use  of  the  applicable  metric  tensor  (Ref. 
10).  Specifically,  Euclidean  3-space  covered  by  the  oblique  Cartesian  coordinate  system 
defined  by  the  transmitted  beam  pattern  of  a multiaxis  LDV  system  is  characterized  by 
a metric  tensor  with  constant  elements.  Therefore,  it  is  reasonable  to  seek  mathematical 
relationships  for  velocity  vector  magnitude  and  spatial  orientation  in  terms  of  the 
applicable  metric  tensor  and  the  LDV-measured  values  of  the  velocity  vector  projections. 

3.0  LDV  MEASUREMENTS  AND  VELOCITY  VECTOR  PROJECTIONS 


Most  present-day  LDV  transmitters  are  of  the  crossed-bean.  configuration  (Ref.  12). 
This  configuration  is  compatible  with  the  reference-beam  receiver  as  well  as  the 
dual-scatter  receiver.  The  expression  for  the  signal  frequency  (Doppler  frequency) 
produced  by  one  axis  of  a crossed-beam,  reference-beam  LDV  system  is 

fD  = (ks  -k,)  (23) 

where  V is  the  velocity  vector  of  a light-scattering  particle  passing  through  the 

— ► 

crossed-beam  or  probe  volume  region,  kj  is  the  propagation  vector  of  the  illuminating 
laser  beam,  and  ks  is  the  propagation  vector  of  the  scattered,  Doppler-shifted  laser 
radiation.  For  nonrelativistic  velocities,  Eq.  (23)  can  be  put  into  the  form 


2 sin  6/2 


(24) 


where  X is  the  wavelength  of  the  illuminating  laser  radiation,  6 is  the  angle  subtended  by 
the  two  intersecting  laser  beams,  and  iff  is  a unit  vector  of  coplanar  with  lc$  and  ltj  and 
perpendicular  to  the  bisector  of  the  angle  0. 
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The  signal  frequency  produced  by  one  axis  of  a crossed-beam,  dual-scatter  type  LDV 
is  given  by 

fD  = 17  V ' ( ki2  -ku)  (25) 

— ► 

where  kn  is  the  propagation  vector  of  the  first  of  two  crossed  beams  of  illuminating  laser 
radiation  and  k|2  is  the  propagation  vector  of  the  second  illuminating  beam.  With  0 as 
the  angle  subtended  by  kt2  and  kI2,  Eq.  (25)  can  be  put  into  a form  identical  to  Eq. 
(24)  which  states  that  an  LDV  system,  whether  reference-beam  or  dual-scatter  type, 
simply  measures  the  orthogonal  projection  of  V upon  the  ut'  axis. 

4.0  THE  METRIC  TENSOR 

As  previously  stated,  a general  coordinate  transformation  is  expressed  by 


= (*‘  ) 


The  total  differential  of  x‘  is  given  by 


dx  = f-  dx 

0x'  (27) 

The  covariant  metric  tensor  is  defined  by  substituting  Eq.  (27)  into  the  expression  for 
scalar  invariance  of  the  infinitesimal  distance  ds  in  rectangular  Cartesian  coordinates  and 
ds'  in  general  curvilinear  coordinates  (Ref.  13): 

(ds')2  = gj'^'dx1  dxj  = (ds)2  = 5- dx1  dxl  (28) 


From  this  it  can  be  seen  that 


0X1  dx> 


For  the  special  case  of  Cartesian  coordinates,  the  general  transformation  given  by  Eq. 
(26)  simplifies  to 


x1  = a:,x’ 

i 


where  A,"'  is  a matrix  of  constants  independent  of  coordinate  position.  Examination  of 
Eqs.  (29)  and  (30)  shows  that  the  elements  of  A|-  must  first  be  determined  before  the 
elements  of  the  metric  tensor  can  be  computed. 
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4.1  COORDINATE  TRANSFORMATIONS  IN  3-SPACE 

In  a Cartesian  coordinate  system  a vector  can  be  represented  as  a directed  line 
segment  connecting  two  points  in  the  space  Such  a vector  is  shown  in  Fig.  2 as  the 
directed  line  segment  OP  and  is  identified  as  V.  The  familiar  parallelogram  law  for  vector 


Figure  2.  Components  of  vector  V with  respect  to  oblique 
Cartesian  coordinate  systems  x'  and  x' . 


addition  of  the  vector  components  with  respect  to  two  different  oblique  Cartesian 
coordinate  systems  is  illustrated  in  Fig.  2.  Further  study  of  this  figure  shows  that  the 
orthogonal  projection  of  V onto  the  x‘  axis  is  given  by 


'1  2 


1 3 


1 1 


X 1 + C j -2  - X 2 


(31) 


where  b,"j"  is  the  cosine  of  the  angle  subtended  by  the  positive  directions  of  the  x1  and 
xi”  axes  and  Cj"j'  is  the  cosine  of  the  angle  subtended  by  the  positive  directions  of  the 
x'"  and  x’’  axes.  In  general. 


(32) 


where  bj'y  and  cpy  are  matrix  arrays  of  these  direction  cosines.  Premultiplication  of 
Eq.  (32)  by  b)"1”  yields 


= bJ 


> j 


(33) 


where  b>'v'  is  the  inverse  of  bj»j".  A similar  procedure  produces 

x>  = b>  * c , y x' 


(34) 
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For  the  special  case  that  xi  = xi,  where  xi  represents  a rectangular  Cartesian 
coordinate  system,  b>"‘"  becomes  the  identity  matrix  Ii‘ 


x'  = E’c.^x'  = A\x' 


which  is  in  agreement  with  Eq.  (30).  Equation  (34)  now  becomes 


x'  = b'  1 c . r x'  = D1  x* 


(35) 


(36a) 


where 


= -J7  jCl'l  E1  (b2'3')  2 ] +c2',  (k1Vb2V-b,v)  +c3'.  (bl'2'b2'3'  ~ b 1 3 ')  J 


(36b) 


Dl  = A'  (b2'3')2]  +C2'2(br3'b2'3'-br2')  +C3'2  (bl'2'b2'3'-‘>1'3')| 


(36c) 


Dl  - X'  jvaf1-  (b2  '3 ' ) 2]  + C2'3  (blVb23'-blv)  +C3'3  (bl'2'b2'3'-b,'3')| 

(36d) 

D1  A7  I C 1 1 (bl'3'b2'3'_bl'2 ')  + tVl  [ 1 - (bl'3')2]  +c3'l  (bl'2'br3'-b2'3')j 


(36e) 


D2  = — 
2 \ ' 


jC.'2(b.Vb2V-b,v)  +C2'2  [*  - (biv)2]  + «V2  (bl'2'b.'3'-b2'3')j 


(360 


D3  - A ' |C»'3  (blVb2'3'-b,  '2')  +C2'3  [l  - (b.'3')2]  +C3'3  (b.  V b.  V " b2 'a')! 


(36g) 


°>  = b |Ci'i(bi'2'b2V-b.v)  (c2'.(b1Vb.Vb2v)*VI  [J-(b.v)2]  ( 


(36h) 
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D!  - r jv, 

(blVb2V-blv)  + C2'2  (bl  '2'bl'3'  "b2'3')  + C3'2 

[>-  (v,)-]| 

' 4-  ]'■'»( 

’br2'b2'3'-  brs')  +C2'3  (b  Tail's'-  b2'3')  + <V3 

(36i) 

.'-(v,-)’]  | 

(36j) 

4-.  1-  (blV) 

I2  ~ (b2V  )2  ~ (bl'3')2  + 2 bl'2/b2VblV 

(36k) 

4.2  ELEMENTS  OF  THE  METRIC  TENSOR 

Expansion  of  Eq.  (35)  yields  the  expressions  for  the  rectangular  coordinates  as  a 
function  of  the  oblique  coordinates.  Differentiation  of  these  expressions  and  substitution 
into  Eq.  (29)  yield  values  for  the  elements  of  the  covariant  metric  tensor: 

*i  V = [(cu'' 

)2+(C-02+(C3r)2] 

(37a) 

82  V = £ (Cl2\ 

)2+  (C22')2+  (C32')2] 

(37b) 

63  V = £(C13'y 

)'+  (^23')  2 + (C33  / ) 2] 

(37c) 

82  '1 ' = 8 1 '2 ' = 

£ (C  11  ' ) (C12')  + (C2l')  ( C 22  / ) + (C3l')  (C32  ')  J 

(37d) 

83  V = 81V  = 

[(  Cir)  (C 13  /)  + (C2l)  (c23  ')  + (c3l')  (C33 ' )] 

(37e) 

83  V = 82  V = 

[(C12')  (C13')  + (C22')  (C23')  + (C32')  (C33') 

| (370) 

5.0  VECTOR  MAGNITUDE 

ExpandinR_JEq.  (14  )and  taking  the  square  root  yield  the  expression  for  the 
magnitude  of  V in  terms  of  the  elements  of  the  covariant  metric  tensor  and  the 
contravariant  measure  numbers: 
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V = [«rr  (vl  )2  + 82'2'(v2  )2  + «3'3'(v3  )2  + 2«1/2'v1  y2 

+ 2«r3'vl  V3  + 2 g2  '3  ' v2  V3  ] ' 


(38) 


Solution  of  Eq.  (17)  for  the  v‘  and  substitution  into  Eq.  (38)  give  the  expression  for 
vector  magnitude  in  terms  of  the  vector  components: 

V = [(V)^(V)^(V)^2glV(glVg2V)-,/1  (V)r(V)2, 

+ 28j  V (biVBsv)  ~/2  (V)r(V)3^  2 g2'3'(82'2'83'3')  “ (V)2'(VV]  '*  (39) 

The  general  expression  for  vector  magnitude  in  terms  of  vector  projections  is  quite 
complicated.  For  practical  applications,  however,  a straightforward  procedure  is  to  solve 
Eq.  (18)  for  the  vi’  in  terms  of  the  (V j_)j'  and  to  then  substitute  the  vi'  into  Eq.  (38). 

.9 

The  vJ  are  given  by 

yl  = \ j(giv)V2(M  rfsaVBsV"  (b2'3  ')  2] 

“ (e2  '2 ')  /2  ( Vi ) 2 ' [«i  '2 ' 83  '3 81  's ' 82  '3 '] 

* ('>  •')  ,/,(Vl)  j (40a| 

y2  = \ j~(glV)  /j  (Vi)r[gi'2'g33'_gi'3'K2'3'] 

+ (82'2'),/2  (\)2'  [8iV83V-(gr3')2  ] 

-(«3v),/2(Vi)3{8lV82V-81V81v]J  (4Qb) 

y3  ' " l j (fix'  r)  HVt).'[8,VB2V  -8,'3'82'2'] 

- (s2v)  /2  (Vl)2'[8l'r82'3'-8r2'8r3'  J 

+ (*3v),/2  (vi) 3' [e ,'1' 82'2'-  (biv)  2 ] J (40c) 


15 


AEDC-TR-79-35 


where 


A = 8rr[g2'2'63'3'-  («2v)  2]  -8l'a'  [bi'2'8s'3'“  8l  '3  ' «2  '3  ' ] 

+ «lV  [biV^VSiv] 


6.0  ANGLE  SUBTENDED  BY  TWO  INTERSECTING  VECTORS 


(40d) 


Dividing  Eq.  (14)  by  V2  yields 

v>  vi 

1 = V ~ “ (41) 

so  that  v*  /V  is  a measure  number  of  the  unit  vector  u?.  Assume  that  a vector  W 
intersects  V to  define  the  angle  0y  ,w  as  shown  in  Fig.  3.  The  unit  vector  uw  and  unit 


0V,W  uvuw 


Figure  3.  Angle  <>v  <w_^ubten^ed  by  the  intersection 
of  vectors  V and  W. 

vector  measure  number  W''/W  are  associated  with  the  vector  W.  The  magnitude  of  the 
vector  (uv  - lfw ) can  be  found  by  the  cosine  law  to  be 


"v  ~H  — COS  ,w) 


This  same  vector  magnitude  is  given  by  [see  Eq.  (24)] 


|-»  -*  | 0 / v1  w1  \ ( v*  w*  \ 

|uv-uw|  = g*'j'yv  W ) \~T~~rJ 


V V'  » V*  _ V w 

Bi'j'  ~v  \T+g'V  ~w  gi V “v  w" 


O-rff) 
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Comparison  of  Eqs.  (42)  and  (43)  shows  that 


^v,w 


arc  cos 


(44) 


6.1  ANGLE  SUBTENDED  BY  V AND  A RECTANGULAR  COORDINATE  AXIS 


The  set  of  measure  numbers  of  a vector  is  a univalent  (rank  one)  tensor  (Ref.  13.)  A 
contra  variant,  univalent  tensor  is  defined  as  a quantity  that  transforms  as  (Ref.  13) 


u 


(45) 


Once  again  let  x‘  represent  a rectangular  Cartesian  coordinate  system,  and  let  x'* 
represent  an  oblique  Cartesian  coordinate  system.  A unit  vector  Uj  in  the  x'  coordinate 
direction  has  a single  non-zero  measure  number  u‘  = 1-  The  corresponding  measure 
numbers  in  the  oblique  coordinate  system  are  given  by 


u 


1,  UJ 


(46) 


Now  Eq.  (44)  can  be  used  to  find  the  angle  subtended  by  the  vector  V and  the  x> 
coordinate  axis: 


<f> 


v,w 


= arc  cos 


/ vi'/\ 

y8>'f  v ) 


(47) 


At  this  point  Eqs.  (36)  and  (45)  can  be  used  to  find  the  u'  in  the  terms  of  u*.  The 
u*  values  can  then  be  substituted  into  Eq.  (47)  to  obtain 

0V  xi  = arcco8jl[glVD;  v1  +g2'2'Df  v2  +83'3'Df  y3 


2'  \ ' 1*2'  3 ' . * 

+ 8l'2'Di  v + <?2'rDi  v + 8 1 '3 ' D i v 


83 


V D ! y3  + 82'3'Di  y2  + «3'2'Df  y3  ] j 


(48) 


The  expressions  for  the  angle  <p  in  terms  of  vector  components  can  be  found  by 
substituting  Eq.  (17)  into  Eq.  (48): 
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<f>y  x2  = arccros  • 


0V,x3  = 


2'  3/  *] 

«1V+D2  *l'2'+D2Sr3'J 

(vy 

(«.v) 1/1 

82V+D2  82'2'  + D2  82v] 

(Vy 

(*2V) 

«3V+D2  «3'2'+D2  »3'3'] 

(V>,' 

(*sv)l/l 

«1V+D3  8l'2'+D3  *1  '3 '] 

<VV 

2‘ ' 3 ' 1 

g2'l'+D3  «2'2'+D3  «2'3'J 

(V>,. 

(«2'2')Vl 

2'  _3 ' 1 

(«3V)/j 

83V+D3  «3'2'+D3  «3 '3 'J 

(49b) 


f (49c) 

The  general,  3-space  expressions  for  0 in  terms  of  the  (VAV  projections  are  quite 
unwieldy.  For  practical  applications,  the  measured  values  of  the  (VAV  can  be  inserted 
into  Eqs.  (40  a through  c);  then  the  resultant  v1'  numerical  values  can  be  substituted  into 
Eq.  (48) to  find  0. 

7.0  COMPONENT  AND  PROJECTION  INTERRELATIONSHIPS 

For  purposes  of  analysis  it  is  often  convenient  to  use  expressions  that  interrelate 
components  and  projections  in  the  x1  and  x1  coordinate  systems.  The  familiar 
parallelogram  law  for  vector  addition  of  components  can  be  used  to  describe  the  vector 
quantity  V in  the  x1  as  well  as  the  x*  coordinate  systems  (see  Fig.  2): 


V = U (V).Z 

i=l  ' ' 


(50) 
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V = I (V)-v  (51) 

i = l 

where  m and  uj'  are  unit  vectors  collinear  with  the  x‘  and  x1'  axes  respectively. 

By  forming  the  scalar  product  of  V as  expressed  in  Eq.  (50)  with  Uj,  one  finds  that 


1 0 0 


(Va)2  =0  1 0 


0 0 1 


which  simply  formalizes  the  earlier  statement  that  (Vx)j  and  (V)j  are  identical  quantities 
in  the  rectangular  Cartesian  coordinate  system  x>.  Forming  the  scalar  product  of  V as 
expressed  in  Eq.  (51)  with  Uj'  yields  the  relationships  for  (V^V  in  terms  of  (V)j' : 


(Vl),'  1 1*  1 '2 ' **  1 '3 ' (^'1' 

(Vi)2'  = b2'l'  1 b2  '3 ' (V)2. 

.(Vi)3'J  Lb3'r  b3'2'  1 J [(V)3'J  (53) 

The  x*  components  in  terms  of  x>’  components  can  be  found  by  forming  the  scalar 
product  of  Eq.  (50)  and  Eq.  (51)  with  U;  and  equating  the  results: 

(V)i  cn'  ci2'  c 13 ' (VV 

(V) 2 = c21'  c22'  c23'  (V) 2 ' 

C31 ' c32'  c33 ' *VVj  (54) 


The  functional  interrelationship  of  (Vi  V and  (VV  can  be  found  by  forming  the 
scalar  product  of  Eq.  (50)  and  Eq.  (51)  with  Uj'  and  equating  the  results: 
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8.0  2-SPACE  RELATIONSHIPS 

The  2-space  versions  of  the  preceding  relationships  are  important  because  of  the 
widespread  use  of  two-axis  LDV  systems  for  two-dimensional  flow  studies.  For 
simplification  let  x3  be  collinear  with  x3  , and  let  both  be  perpendicular  to  x1  , x2\x! 
and  x2.  Then,  for  2-space  we  have 

v3  = (V)3'  = (Vi),.  = c31'  = c32'  = = c3'2 


= D1  = D2'=  D3' 
3 3 1 

= °2  = °3  = g i '3 ' = «2'3'  = 0 

(57) 

It  can  also  be  readily  shown  that 

c3'3  = c33 ' = 

g2'2'=  1 

(58) 

gj'2' 

- b j ^ ' 

(59) 

D1  “ .in 2y  [ 

cn  ~ c2'i  bi'2^] 

(60) 
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D2 


A 

•2,, 

[C1'2  - c2'2  bl'2'] 

. - 

sin  y 

1 

. 2 

sin  y 

- cri  b j '2 ' + c 2 'j 

1 

, ”1 

sin2y  | 

“ C1'2  bl'2'+  c2'2  J 

(61) 


(62) 


where  7 = arccos  biV. 

8.1  VECTOR  MAGNITUDE 

In  2-space,  Eq.  (38)  becomes 

v-  [(«*'K  (»2-)2*2b,v,'»2'],/- 

while  the  expression  for  V in  terms  of  vector  components,  Eq.  (39),  reduces  to 
V = [(V)  j,  + (V)  ^ + 2blV(V)r(V)2,J 


(63) 

(64) 


(65) 


r 2 ' ' 1 


(66) 


Substitution  of  Eqs.  (58)  and  (59)  into  Eqs.  (40a),  (40b),  and  (40d)  yields 

[('  i)  1"  b>'2‘  (Vi)  2'] 

r=^[-b.-’-(V.)^  (M*'] 


(67) 


(68) 


Now  Eqs.  (67)  and  (68)  can  be  substituted  into  Eq.  (65)  to  find  the  simplified  2-space 
expression  for  V in  terms  of  vector  projections: 


).i n2y  [ (Vl)  l'+  (Vi)  2'~2bl'2'(V1)  1'  (VJ  2']  | 2 


(69) 


It  can  be  seen  that  Eqs.  (66)  and  (69)  both  reduce  to  the  familiar  form  for 
rectangular  Cartesian  coordinates  when  7 = 90  deg: 


v = | (v)j.  ♦ (v 


(70) 


J 
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8.2  THE  ANGLE  <t> 


In  2-space,  the  simplified  expressions  for  <p  in  terms  of  vector  components  are  found 


<AV,xl  = [«Vl  (V)l'+  c2'l  <V)2']} 

0V>X2  = arccosj^  [c  r2  vV)r  + (V>2]| 


The  expressions  for  <p  in  terms  of  velocity  projections  are  given  by 


<£y  xi  = arccos 


\~[  (cl'l  -c2'lbl'2')  (Vj)r  + (c2'l-cribl'2')(Vl)2/]J 


^V,x2  = arccos)  — I (cj'2  - c2-2  bi'2')  (VA)  i ' + (c2'2  “ C1'2  bl'2')  C Vx)  2 ' | f 
(Vainly  J " 


It  can  be  seen  that  Eqs.  (71)  and  (73)  and  Eqs.  (72)  and  (74)  reduce  to  the  familiar 
forms  for  rectangular  Cartesian  coordinates  when  y = 90  deg: 

<£V>X1  = arccos  | i j^cri  (,V)r  + cn(V)2'  ^ 

and 

<£v,x!  = “rccos  j^[cj'2(v)r  + c2.2(V)2.J  | ^ 

8.3  COMPONENTS  A»U  PROJECTIONS 

In  2-space,  Eqs.  (52)  through  (56)  reduce  to 


1 0 


0 1 
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9.0  CONCLUSIONS 

The  complexity  of  the  general  3-space  relationships  clearly  illustrates  the  desirability 
of  implementing  orthogonal  LDV  measurement  axes  collinear  with  the  reference  system 
axes.  When  these  conditions  are  unattainable,  the  equations  in  this  report  can  be  used  to 
compute  the  values  for  velocity  magnitude  and  spatial  orientation  within  the  uncertainties 
dictated  by  the  remaining  LDV  system  error  mechanisms. 

The  provision  of  the  correct  equations  to  account  for  axis  nonorthogonality  and 
rotation  now  focuses  attention  upon  the  status,  with  regard  to  accuracy  and  precision,  of 
techniques  and  devices  for  determining  the  bjy  and  cfj  values  to  be  used  in  these 
equations.  Moreover,  further  reductions  in  LDV  instrument  errors  will  require  additional 
refinement  of  (1)  techniques  and  devices  for  determining  the  signal  frequency-to-velocity 
calibration  factor  for  each  LDV  measurement  axis  and  (2)  techniques  and  devices  for 
measuring  the  signal  frequency  associated  with  each  axis. 

REFERENCES 

1.  Rabinawitz,  P.,  Jacobs,  S.,  Targ,  R.,  and  Gould,  G.  “Homodyne  Detection  of 
Phase-Modulated  Light.”  Proceedings  of  the  IRE,  Vol.  50,  No.  11,  November  1962, 
p.  2365. 


23 


AEOC-TR  79  35 

2.  Cummins,  H.  Z.  and  Knable,  N.  “Single  Sideband  Modulation  of  Coherent  Light  by 
Bragg  Reflection  from  Acoustical  Waves.”  Proceedings  of  the  IEEE,  Vol.  51,  No.  9, 
September  1963,  p.  1246. 

3.  Cummins,  H.Z.,  Knable,  N.,  and  Yeh,  Y.  “Observation  of  Diffusion  Broadening 
Rayleigh  Scattered  Light.”  Physical  Review  Letters,  Vol.  12,  No.  6,  10  February 

1964,  pp.  150-153. 

4.  Yeh,  Y.  and  Cummins,  H.  Z.  “Localized  Fluid  Flow  Measurements  with  an  He-Ne 
Laser  Spectrometer.”  Applied  Physics  Letters,  Vol.  4,  No.  10,  15  May  1964,  pp. 
176-178. 

5.  The  Accuracy  of  Flow  Measurements  by  Laser  Doppler  Methods.  Proceedings  of  the 
LDA-Symposium.  Copenhagen,  Denmark,  August  25-28,  1975. 

6.  Crosswy,  F.  L.  and  Brayton,  D.  B.  “Design  and  Evaluation  of  a Laser  Doppler 
Velocimeter  for  the  AEDC-PWT  4-ft  Transonic  Tunnel.”  AEDC-TR-78-30,  July 
1978. 

7.  Rolfe,  E.,  Silk,  J.  K.,  Booth,  S.,  Meister,  K.,  and  Young,  R.M.  “Laser  Doppler 
Velocity  Instrument.”  NASA  CR-1199,  December  1968. 

8.  Barnett,  D.  O.,  and  Giel,  T.  V.  “Application  of  a Two-Component  Bragg-Diffracted 
Laser  Velocimeter  to  Turbulence  Measurements  in  a Subsonic  Jet.”  AEDC-TR-76-36 
(AD-A025355),  May  1976. 

9.  Stratton,  Julius  Adams.  Electromagnetic  Theory,  McGraw-Hill  Book  Company,  Inc., 
New  York,  1941 . 

10.  Sokolnikoff,  Ivan  Stephen.  Tensor  Analysis:  Theory  and  Applications  to  Geometry 
and  Mechanics  of  Continua.  John  Wiley  and  Sons,  Inc.,  New  York,  1964. 

11.  Lass,  Harry.  Vector  and  Tensor  Analysis.  McGraw-Hill,  New  York,  1950. 

12.  Lennert,  A.  E.,  Brayton,  D.  B.,  and  Crosswy,  F.  L.  “Summary  Report  of  the 
Development  of  a Laser  Velocimeter  to  be  Used  in  AEDC  Wind  Tunnels.” 
AEDC-TR-70-101  (AD871321),  July  1970. 

13.  Moon,  Parry  and  Spencer,  Domina  Eberle.  Vectors.  Van  Nostrand,  Princeton,  N.  J., 

1965. 


24 

i 

1^ 


AEDC  TR  79  35 


A'.' 


ap 


bj'i' 


b)'1' 


ci'j 

-*■ 

dr 

ds 

ds’ 


NOMENCLATURE 

Matrix  array  of  constants,  Eq.  (30) 

Unitary  vector,  Eq.  (4) 

Reciprocal  unitary  vector,  Eq.  (6) 

Matrix  array  of  direction  cosines,  Eq.  (34) 

Inverse  of  bjV,  Eq-  (36a) 

Matrix  array  of  direction  cosines,  Eq.  (36a) 

Matrix  array  of  constants,  Eq.  (36a) 

Differential  displacement  vector,  Eq.  (3) 

Differential  displacement'  magnitude  in  the  x‘  coordinate  system,  Eq. 
Differential  displacement  magnitude  in  the  x*'  coordinate  system,  Eq. 
I T)V  siunal  frenneiHW.  En.  03 1 


(22) 

(12) 


gi'j' 

.t.t 

g'J 

k* 

-► 

ki 

— ► 

r 

3 


uv 

-► 

uw 


c 

V 


Covariant  metric  tensor,  Eq.  (12) 

Contravariant  metric  tensor,  Eq.  (13) 

Scatter  radiation  propagation  vector,  Eq.  (23) 

Illuminating  beam  propagation  vector,  Eqs.  (23)  and  (25) 
Position  vector,  Eq.  (3) 

Unit  vector  collinear  with  the  x‘  axis,  Eq.  (50) 

Unit  vector  collinear  with  the  xj'  axis,  Eq.  (24) 


